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separate coupled cavities via transitionless quantum driving 

Wu-Jiang Shan 1 , Ye-Hong Chen 1 , Yan Xia 1 ’*, and Jie Song 2 ’* 

1 Department of Physics, Fuzhou University, Fuzhou 350002, China 
2 Department of Physics, Harbin Institute of Technology, Harbin 150001, China 

By jointly using quantum Zeno dynamics and the approach of “transitionless quan¬ 
tum driving (TQD)” proposed by Berry to construct shortcuts to adiabatic passage 
(STAP), we propose an efficient scheme to fast generate multiatom Greenberger- 
Horne-Zeilinger (GHZ) states in separate cavities connected by opitical fibers only 
by one-step manipulation. We first detail the generation of the three-atom GHZ 
states via TQD, then, we compare the proposed TQD scheme with the traditional 
ones with adiabatic passage. At last, the influence of various decoherence factors, 
such as spontaneous emission, cavity decay and fiber photon leakage, is discussed by 
numerical simulations. All the results show that the present TQD scheme is fast and 
insensitive to atomic spontaneous emission and fiber photon leakage. Furthermore, 
the scheme can be directly generalized to realize iV-atom GHZ states generation by 
the same principle in theory. 
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I. INTRODUCTION 


Quantum entanglement is not only one of the most important features in quantum me¬ 


chanics |lt, but also a key resource for testing quantum mechanics against local hidden theory 


Recently, the entangled states have been applied in many fields in quantum informa¬ 
tion processing (QIP), such as quantum computing |^|, quantum cryptography [4], quantum 
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teleportation |s, 6|, quantum secret sharing Q, and so on. These promising applications 
have greatly motivated the researches in the generation of entangled states. 

It is worth noting that a typical entangled state so-called Greenberger-Horne-Zeilinger 
(GHZ) state | GHZ) = -^=(|000) + |111)) , hrst proposed and named by Daniel M. Green- 
berger, Machael Horne and Anton Zeilinger 8|, has raised much interest. Contrary to other 
entangled states, the GHZ state exhibits some special features, such as it is the maximally 
entangled state and can maximally violate the Bell inequalities [ 9 ] . In 2001, Zheng has pro¬ 
posed a scheme to test quantum mechanics against local hidden theory without the Bell’s 


inequalities by use of multiatom GHZ state 


10]. Therefore, great interest has arisen regard¬ 


ing the significant role of GHZ state in the foundations of quantum mechanics measurement 
theory and quantum communication. At present, the hrst and main problem we face is how 
to generate GHZ state by using current technologies. To our knowledge, in some experimen- 


tal systems, such as trapped ions systems U , photons systems HQ, and atoms systems 
[141 ]. scientists have realized the generation of such GHZ state. Recently, a promising ex¬ 
perimental instrument named cavity quantum electrodynamics (C-QED), which concerns 
the interaction between the atom and the quantized held within cavity 15j], has aroused 
much attention. Based on C-QED, many theoretical schemes for generating GHZ state have 
been proposed. For example, Li et al. have proposed a scheme to generate multiatom GHZ 


state under the resonant condition by Zeno dynamics 


16], but the scheme is sensitive to 


the atomic spontaneous emission and hber photonic leakage. Hao et al. have proposed an 
efficient scheme to generate mulitiatom GHZ state under the resonant condition via adia¬ 
batic passage [17j, but it takes too long time. Chen et al. have proposed a smart scheme to 
overcome the above drawbacks, but the scheme needs to trap three atoms in one cavity 18], 
such design is difficult to manipulate each atom in experiment and to construct a large-scale 
quantum network. 

On the other hand, in modern quantum application held, an important method to ma¬ 
nipulate the states of a quantum system is adiabatic passage, included “rapid” adiabatic 
passage (RAP), stimulated Raman adiabatic passage (STIRAP), and their variants 19]. 
The adiabatic passage covers the shortage with respect to errors or fluctuations of the pa¬ 
rameters compared with the resonant pulses, but its evolution speed is very slow, so it may 
be useless in some cases. In recent years, shortcuts to adiabatic passage (STAP), which 
accelerates a slow adiabatic quantum process via a non-adiabatic route, has aroused a great 
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deal of attention. Many 


population transfer 


heoretical proposals have been presented to realize QIP, such as fast 


quantum phase gates 


20-2 


23L fast entanglement generation 22, 24], fast implementation of 


26], and so on. To our knowledge, the main methods to construct 


effective shortcuts has two forms: one is invariant-based inverse engineering based Lewis- 


Riesenfeld invariant (IBLR) |27j and the other is transitionless quantum driving (TQD) 


which is pointed out by Berry. The two methods are strongly related 


29], but also have 


their own characteristics. For example, the former does not need to modify the original 
Hamiltonian H 0 (t), but the algorithm is suitable for some special physical models. The lat¬ 
ter needs to modify the original Hamiltonian H 0 (t) to the “counter-diabatic driving” (CDD) 
Hamiltonian H(t) to speed up the quantum process. The fixed Hamiltonian H(t) can be 
obtained in theory, but it does not usually exist in real experiment. 

In addition, the quantum Zeno effect (QZE), first understood by Neumann 30] and 
named by Misra and Sudarshan 31], exhibits a especially experimental phenomenon that 
transitions between quantum states can be hindered by frequent measurement. The system 


will evolve away from its initial state and remain in the so-called “Zeno subspace” c 


the measure due to frequently projecting onto a multi-dimensional subspace 32 


ehned by 


33]. This 


is so-called quantum Zeno dynamics (QZD). Without making using of projection operators 
and non-unitary, “a continuous coupling” can obtain the same quantum Zeno effect instead 


of discontinuous measurements 


M 


35]. Now, we give a brief introduction of the quantum 


Zeno dynamics in the form of continuous coupling 35]. Suppose that the system and its 
continuously coupling external system are governed by the total Hamiltonian H to t = H s + 
KH ei where H s is the Hamiltonian of the quantum system to be investigated, H e is an 
additional Hamiltonian caused by the interaction with the external system, K is the coupling 
constant. In the limit K —> oo, the evolution operator of system can be expressed as 


U(t) = exp [-it Y,n( Kr lnPn + P n H a P n )\, W 


rere P n is the eigenprojection of H e corresponding 


to the eigenvalue J] n , i.e., H e P n = rj n P n 36]. 


Inspired by the above useful works, we make use of Zeno dynamics and TQD to construct 
STAP to generate iV-atom GHZ state in C-QED. Our scheme has the following advantages: 
(1) The atoms are trapped in different cavities so that the single qubit manipulation is more 
available in experiment. (2) The fast quantum entangled state generation for multiparticle 
in spatially separated atoms can be achieved in one step. (3) Numerical results show that our 
scheme is not only fast, but also robust against variations in the experimental parameters 




































4 


and decoherence caused by atomic spontaneous emission and fiber photon leakage. In fact, 
further research shows that, the total operation time for the scheme is irrelevant to the 
number N of qubits. 

The paper is organized as follows. In section II, we give a brief introduction to the 
approach of TQD proposed by Berry. In section III, we introduce the physical modal and 
the systematic approximation by QZD. In section IV, we propose the schemes to generate 
the three-atom GHZ state via TQD and adiabatic passage, respectively. The decoherence 
caused by various factors is discussed by the numerical simulation. In section V, we directly 
generalize the scheme in section IV to generate iV-atom GHZ states in one step. At last, we 
discuss the experimental feasibility and make a conclusion about the scheme in section VI. 

II. TRANSITIONLESS QUANTUM DRIVING 

Suppose a system is dominated by a time-dependent Hamiltonian H 0 (t ) with instanta¬ 
neous eigenvectors \4> n (t )) and eigenvalues E n (t ), 

H 0 (t)\Mt)) = E n (t)\Mt))- (1) 

When a slow change satisfying the adiabatic condition does, the system governed by H 0 {t) 
can be expressed at time t 

1 = e%Ut) \ 0n(*))> 

£n(t) = J o dt'E n (t')+i J dlf(<f> n (lf)\dt><j) n (lf)), ( 2 ) 

where d t > — Because the instantaneous eigenstates | <j) n (t)) do not meet the Schrodinger 
equation ihd t \4> n {t)) = H 0 \(j> n (t)), a finite probability that the system is in the state |d m / n (f)) 
will occur during the whole evolution process even under the adiabatic condition. 

To construct the Hamiltonian H(t) that drives the instantaneous eigenvector 1 4> n (t)) ex¬ 
actly, i.e., there are no transitions between different eigenvectors during the whole evolution 
process, we define the unitary operator 

U = (3) 

n 

we can formally solve the Schrodinger equation 


H(t) =ih(d t U)Ul 


( 4 ) 
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Substituting eq. (3) into eq.(4), the Hamiltonian H(t) can be expressed 

H{t) = {\d t <f) n )((t) n \ - 11^2 |0n)6»(0n|), (5) 

n n 

the simplest choice is E n = 0, for which the bare states 1 (j> n (t)), with no phase factors, are 
driven by 

H(t) = I d t (f>n){4>n\- (6) 


III. PHYSICAL MODAL AND SYSTEMATIC APPROXIMATION BY QZD 


For the sake of the clearness, let us first consider the physical modal that three identical 
atoms ai, o >2 and a 3 are trapped in three linearly arranged optical cavities C \, C 2 and C 3 , 
respectively. As shown in FIG. [I], each atom possesses one excited level |e) and three ground 
states | gi), \g 0 ) and \g r ). The cavities C\ and C 3 are single-mode, the cavity C 2 are bi-mode. 
Ci, C 2 and C 3 are connected by the optical fibers / 1 , / 2 , respectively. Assuming that the 
transition |e) ai(3) O | g 0 )a 1(3) is resonantly driven by a external classical field with the time- 
dependent Rabi frequencie 0 1 ( 3 )()f:), while the transition |e)i( 2 ) ^ |^?z) 1 ( 2 )(|e) 2 ( 3 ) -H- | 5 ' r ) 2 ( 3 )) 
is resonantly coupled to the left-circularly(right-circularly) polarized cavity mode with the 
coupling constant gi( r ), respectively. 

In the short-fiber limit, i.e., (2Lh)/(27rc) 1 (L is the length of the fibers, v is the decay 

rate of the cavity fields into a continuum of fiber modes and c is the speed of light), only 
one resonant mode of the fiber interacts with the cavity mode 37]. In the rotating frame, 
the Hamiltonian of the whole system can be written as (U = 1) 


Htot — Hi + H c , 

Hi = ^ n o (t)\e) ao (go\ +H.c., 

0=1,3 

He = gildll\e) ai (gi\ + 92l^2l\&) a 2 {9l\ + 92r&2r\e) a 2 {9r\ + 93r^3r\^) a 3 {9r 


+vib\(au + a 2 i) + V 2 blia 2r + a 3r ) + H.c., 


( 7 ) 
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where a' kl (a£ r ) and aki {a>kr) denote the creation and annihilation operators for the left- 
circularly (right-circularly) polarized mode of cavities Ck (k = 1, 2, 3), respectively; bj and 
bj denote the creation and annihilation operators associated with the resonant mode of fiber 
fj (j = 1, 2), respectively. For the sake of simplicity, we assume gu = g 2 i = g- 2 r = g%r = g 
and V\ = v 2 = v. If the initial state of the whole system is \g 0 gig r ) |0)ci |00)c 2 |0)c 3 |0)/i |0)/ 2 
(here | g 0 9i9r) — \9o9i9r)aia 2 a 3 )i the whole system evolves in the following subspaces 

I0i) = l5 , o5 , i5 , r)|0)ci|00)c 2 |0)c 3 |0)/ 1 |0) /2 , |0 2 ) = |e5 , ;fi'r)|0)ci|00)c 2 |0)c 3 |0)/ 1 |0) /2 , 

103) = lfl , ;fl , i5 , r)|l)ci|00)c 2 |0)c 3 |0)/ 1 |0)/ 2 , |0 4 ) = l5 , ifl , i5'r)|0)c 1 |00)c 2 |0)c 3 |l)/i|0)/ 2 , 

105) = lfl , «fl , i5 , r)|0)ci|10)c 2 |0)c 3 |0)/ 1 |0)/ 2 , 10 6 ) = |5 , ;e(yf T .)|0)c 1 |00)c 2 |0)c 3 |0)/ 1 |0)/ 2 , 

107) = lfl , «5 , r5 , r)|0)ci|01)c 2 |0)c 3 |0)/ 1 |0)/ 2 , |0 8 ) = \gWr9r) |0)ci |00)c 2 |0)c 3 10)1l)/ 2 , 

|09> = lfl , ;5 , r5 , r)|0)c 1 |00)c 2 |l)c 3 |0)/ 1 |0)/ 2 , 1 0io) = |^re)|0)c 1 |00)c 2 |0)c 3 |0)/ 1 |0)/ 2 , 

10ii) = lfl , ;5 , r5 , o)|0)ci|00)c 2 |0)c 3 |0)/ 1 |0) /2 , (8) 

where \ijk) (i, j, k G [e, gi, g Q , g r ]) denotes the state of the atoms in every cavity, |n) s (s = 
Ci, C' 3 , /i, f 2 ) means that the quantum held state of system contains n photons. | n\n 2 )c 2 
means that the number of left-circularly photon is n\ and the number of right-circularly 
photon is n 2 in the cavity C 2 . 

Under the Zeno condition g,v S> , U2 3 , the Hilbert subspace is split into nine invariant 
Zeno subspace 

Zi = {I0l)> |0l), |01l)}) Z 2 = {|0 2 )}, 

Z 3 = (|0 3 )}, ^4 = {104)}, ^5 = {|05>}, Zq = {|0g)}, 

Z 7 = (|0 7 )>, Z 8 = { |0 8 )}, Z 9 = { |0 9 )}, 
where the eigenstates of H c are 
\*Pl) = JVl(|02> - + 106) - ||<M + I'M), 


(9) 
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102) — ^2 (— 102) + £l|</>3) — 'hi 104) _ Xl 105) + Xl 107) + '^l 108) — e l|09) + |0io))j 

103) = ^V 3 (—10 2 ) - £i| 03) - hi|04) + Xi|05> - Xi|0r) + ?7i|0s> + ei|0 9 ) + |0io)), 

1 ^ 4 ) = ^ 4 ( 1 ^ 2 ) — A *1 1^ 3 ) — Cl 104) + ^1 105) — ^ 1 106) + ^1 10r) — Cl 108) — A *1 1 ^ 9 ) + |0io)) 1 

|05) — -^ 5 ( 102 ) + A *1 10 3 ) — Cl 104) — ^1 105) — ^1 106) — ^1 107) — Cl 108) + /^1 109) + |0io))j 

|0e) = ^V 6 ( —102> + S2|03> - 72 104) + X205) - X2|07) + 72 108) “ C 2 |0 9 ) + |01O>), 

|0 7 ) = iV 7 (—|0 2 ) — ^2 10 3 ) — '72104) — X 2 I 05 ) + X 2 I 07 ) + 72 108) + ^2 109) + |0io))j 

|08) = ^Vg (1 02 ) — A* 2 10 3 ) + C 2 104) — &2 | 05) + ^ 2 106) — ^2 |07) + C 208 ) — /^-2 10 9 ) + |0io))j 

|09) = ^Vg(102) + A*210 3 ) + C 2 104) + <^2105) + ^2 106) + <^2107) + 02108) + A*210 9 ) + |01o)X0O) 

with the corresponding eigenvalues 

Ai = 0, A 2 = —\J(g 2 + 2v 2 — A)/2, A 3 = vV + 2u 2 - A)/2, 

A 4 = — a/ (3g 2 + 2v 2 — A)/2, A 5 = \/(3 g 2 + 2v 2 - A)/2, X 6 = -^(g 2 + 2v 2 + A)/2, 


A 7 — \J ( g 2 + 2u 2 + A)/2, A§ — — a/(3 g 2 + 2v 2 + A)/2, Ag — a/(3 g 2 + 2v 2 + A)/ 2, (11) 


where the parameters are 


J g 2 + 2v 2 — A —g 2 + 2v 2 — A J g 2 + 2v 2 — A(g 2 + A) 

e i = - t= -, m = -—-, Xi = 


V2g 


2 gv 


2 V2gi 


\j3g 2 + 2v 2 — A —g 2 — 2v 2 + A a /3 g 2 + 2v 2 — A(—g 2 + A) —g 2 + A 

hi = - 7^ -:> Cl — -iT - -> m — - 7T~r ^—n-> m = 


V2g 


2 gv 


2\f2gi 


J g 2 + 2v 2 + A —g 2 + 2v 2 + A J g 2 + 2v 2 + A(—g 2 + A) 

£2 = - t= -, V2 =-v-, X2 =- 


V2g 


2 gv 


2V2gi 


\/3g 2 + 2v 2 + A g 2 + 2v 2 + A s/3g 2 T2v 2 ~+A(g 2 + A) g 2 + A 

h 2 = - k -> C 2 =--, S 2 =- „ ^ 0 -, 0 2 = (12) 


V2g 


2 gv 


2\[2gi 
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in addition, A = \Jg A + 4n 4 and N w is the normalization factor of the eigenstate \ip w ) (w = 
1,2,•••,9). 

The projector in the kth Zeno subspace Z k is 

pH = \mi (1/3) e z„). ( 13 ) 

The Hamiltonian in Eq. (8) can be approximately given by 

H„ ~ YH >*PH + P^HiPH 

/c,/3,7 

9 

= Afc|V>fe)(V>fc| + Ni(Qi\(f>i)(^i\ + O 3 |0n)(^i| + H.C.). (14) 

k =2 

If the initial state is |5 , o fi , «5 , r)|0)c'i|00)c 2 |0)c 3 |0)/ 1 |0)/ 2 , it reduces to 

H e ff = ^i(^i|0i)(t/ ; i| + ^ 3 10H) (f/h I + H.c. ), (15) 

which can be treated as a simple three-level system with an excited state |-0i) and two 
ground states \4>i) and |0n). Then we obtain the eigenvectors and eigenvalues of the effective 
Hamiltonian H e ff as 


I Vo(t)) = 


cos 9{t) 


— sin 9(t) 


, \v±(t)) = 


sin 9(t) 


cos 9{t) 


with the corresponding eigenvalues rj o = 0 and rj± = and tan 9 


^ and Q = 


IV. THE GENERATION OF THE THREE-ATOM GHZ STATE VIA 
TRANSITIONLESS QUANTUM DRIVING AND ADIABATIC PASSAGE 

A. Adiabatic passage method 

For the sake of the clearness, we first briefly present how to generate the three-atom GHZ 
state via adiabatic passage. When the adiabatic condition \(n 0 \d t n±)\ <C |A' ± | is fulfilled well 
and the initial state is |^(0)) = \<j>i), the state evolution will always follow |n 0 (t)) closely. 
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To generate the three-atom GHZ states via the adiabatic passage and meet the boundary 
conditions of the fractional stimulated Raman adiabatic passage (STIRAP), 


lim 

£—>■—oo 


Hj(t) 

n 3 (t) 


0, lim 


n 3 (t) 


t —»+oo (t) 


tan Ol 5 


(17) 


we need properly to tailor the Rabi frequencies fii(t) and H 3 (t) in the original Hamiltonian 

Htot 


ih{t) = sin a \!q exp [- — -J, 

K 

n -\ o f — (t + t 0 — tf/2) 2 f — (t — t 0 — tf/2) 2 

where Ho is the pulse amplitude and tf is the operation time. t c and t 0 are some related 
parameters to be chosen for the best performance of the adiabatic passage process. In 
order to achieve better performance and meet the boundary conditions, we suitably chose 
the parameters that tana = 1, t Q = 0.14 tf and t c = 0.19 tf. As shown in Fig. 2, the 
time-dependent Hi(t)/Ho and H 3 (t)/H 0 versus t/tf are plotted with the fixed values t 0 and 
t c . With the above parameters, we obtain our wanted three-atom GHZ state ^(t/)) = 
(|0i) — |0n))/\/2 via the adiabatic passage. But this evolution process needs a relatively 
long time to satisfy the adiabatic condition. We will detail the reasons in the section of 
numerical simulations and analyses. 


B. Transitionless quantum driving method 


To reduce the evolution time and obtain the same state as the adiabatic passage, we 
use the approach of TQD to construct STAP. As introduced in the above, STAP speeds 
up a slow adiabatic passage via a non-adiabatic passage route to achieve a same outcome, 
and the TQD method is a important route to construct shortcuts. According to the ideas 
proposed by Berry 28], the instantaneous states in Eq. (16) do not meet the Schrodinger 
equation, i.e., id t \rik) ^ H e ff\nk)(k = 0 ,±), so the situation that the system starts from 
the state |^ n (0)) and ends up in the state |^ m ^ n (t)) occurs in a finite probability even 
under the adiabatic condition. To drive the instantaneous states | nk)(k = 0, ±) exactly, we 
look for a Hamiltonian H(t) related to the original Hamiltonian H e ff according to Berry’s 
transitionless tracking algorithm 28]. From section II, we know the simplest Hamiltonian 
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77(f) possessed the form, 

H(t) = \d t n k (t))(n k (t)\. (19) 

o,± 

Substituting Eq. (16) in Eq. (19), we obtain 

H(t)=ie\<f> 1 )(<f> 11 \+H.c.,. (20) 

where 9 = [Oi (t )03 (t) — Oi (t)Os(t)]/ri 2 . This is our wanted CCD Hamiltonian to construct 
STAP, and we will detail how to construct this Hamiltonian in experiment later. 

For the present system, the CDD Hamiltonian 77(f) is given in Eq. (18), but it is 
irrealizable under current experimental condition. Inspired by Refs. [17, 19], we find an 
alternative physically feasible (APF) Hamiltonian whose effect is equivalent to Hi(t). The 
design is shown in Fig. 3, the atomic transitions is not resonantly coupled to the classical 
lasers and cavity modes with the detuning A. The Hamiltonian of the system reads H' tot = 
H c + Hi + Hd, where Hd = Y^t=\ A|e)fc(e|. Then, similar to the approximation by QZD in 
section III, we also obtain an effective Hamiltonian for the non-resonant system 

K ff = I +n 3 |0n)(^i| + H.C.) +3A7V 1 2 |V> 1 )(Vh|. (21) 

When the large detuning condition 3A Ni > TR, fl 3 is satisfied, we can adiabatically elimi¬ 
nate the state |^i) and obtain the final effective Hamiltonian 

Hfe = ” 7^-101) (011 - 2 ^| 01 l)( 01 l|-^^(I0l)(01l| + 101l) (011) • (22) 

For simplicity, we set = C2 3 = 72(f). The front two terms caused by Stark shift can be 
removed and the Hamiltonian becomes 


77(f) = O00 1 )(0 11 |+tf.c., 


(23) 


where Q x = — 7 ^. The equation has a similar form with Eq. (20), but the effective couplings 
between i,9 and 72^ exist 37 r/ 2 -dephased. To guarantee their consistency, we put a change 
that 72 3 —> —i72 3 . Then, the eigenstates of H e ff become 


Wait)) 


^ cos 6{t) ^ 
0 

^isind(f) J 


, W±(t)) 


1 

71 


^ sind(f) ^ 

±1 

^ — i cos 6[t) j 


(24) 
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and the CDD Hamiltonian H{t) becomes 


H(t) = - 0\4> u )(<h 


(25) 


Compared Eq. (23) with Eq. (25), we can easily get the CDD Hamiltonian when the 
condition D x = — 6 is satisfied. 


fil(i) = a : 3 (() = n(i) = vT/V). 


(26) 


C. Numerical simulations and analyses 

Next we will show that it takes less time to get the target state on the situation governed 
by the APF Hamiltonian H' tot via TQD than by the original Hamiltonian H to t via adiabatic 
passage. The time-dependent population for any state \ij}) is defined as P = |('0|p(f)|'0}|, 
where p{t) is the corresponding time-dependent density operator. We present the fidelity 
versus the laser pulses amplitude D 0 and the operation time t/tf via adiabatic passage in 
Fig. 4. As shown in Fig. 4, we can know that the bigger the laser pulse amplitude is, 
the less time that the system evolution to the target state needs. However, we need to 
satisfy the Zeno conditions j,d> Di,D 3 , so we set Do = 0.2 g. In Fig. 5, we display the 
time-dependent populations of the states \4>i), It/’ target ), and |0n) via adiabatic passage. As 
depicted in Fig. 4 and Figs. 5, the operation time needs tf > 400 /g to achieve an ideal 
result at least. It is awkward in some case. 

Next we will detail the evolution governed by the APF Hamiltonian H' tot via TQD. Ac¬ 
cording to eq. (24) we finally get a GHZ state |= ^(|0i) + *|</>n))- I n Fig. 6, we 
present the relationship between the fidelity of the three-atom GHZ state (governed by the 
APF Hamiltonian) and two parameters A and tf when D 0 = 0.2 g to satisfy the Zeno condi¬ 
tion, where the fidelity of the three-atom GHZ state is defined as F — \(GH Z\p{t f)\GH Z)\ 
(p(tf) is the desity operator of the whole system when t = tf). We find that a wide range 
for parameters A and tf can obtain a high fidelity of the three-atom GHZ state, and the 
fidelity increases with the increasing of A and the decreasing of tf. In order to satisfy the 
large detuning condition, we set A = 2.3 g. The Fig. 6 reveals that the operation time needs 
tf > 72 /g via TQD at least. In Figs. 7 we plot the operation time for the creation of the 
GHZ state governed by H' tot and by H tot with the parameters that tf = 72/g, D 0 = 0.2 g, 
A = 2.3 g and g = v. Numerical results show that the APF Hamiltonian H' tot can govern 
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the evolution to a perfect GHZ state \^{tf)) from l^) in a relatively short interaction time 
while the original Hamiltonian H tot can not. 

In above analysis, we do not consider the influence of decoherence caused by various 
factors, such as spontaneous emissions, cavity decays and fiber photon leakages. In fact, the 
decoherence is unavoidable during the evolution of the whole system in experiment. The 
master equation of the whole system is written as 

P = -i[H tot ,p\ 

3 

+ y( 2 cr fc P°t ~ a t a kP - P°t°k) 

k =1 
2 

+ - a tk a i,kp - P a t,k a i,k) 

k =1 
3 

+ 'y ^ ——(2a 7 ',kP^ k — — P a \ ) k 0j i'yk) 

k= 2 
2 

+ ^2^y( 2b kpK ~ b kb k p- pK b k), ( 27 ) 

k =1 

where yy is the atomic spontaneous emission rate for the fctli atom and K c (f) is the decay 
rate of the kth cavity (/cth fiber), oy denotes the atomic transition from the ground states 
| m) (m = g 0 , gi , g r ) to the excited state |e). For the sake of simplicity, we assume that 
7 i = 72 = 73 = 7 , ki Cl = «c 2 = k C3 = k c and k./, = Kf 2 = Kf. As shown in Fig. 8 , we 
plot the fidelity governed by the APF Hamiltonian H[ ot and by the original Hamiltonian 
H tot and the dimensionless parameters y/g, K c /g and Kf/g, respectively. We can draw a 
conclusion that the fidelities are almost unaffected by the fiber decay both via TQD and via 
adiabatic passage. We focus on the main decoherence factors included the cavity decay and 
the atomic spontaneous emission. As shown in Figs. 9, we plot the fidelity versus the cavity 
decay and the atomic spontaneous emission. We can know the most important decoherence 
factor is the cavity decay. This result can be understood from Ref. [18] that if the Zeno 
condition can not be satisfied very well, the populations of the intermediate states including 
the cavity excited states can not be suppressed ideally. 

From the above anslysis, we can obviously know that the evolution time from the initial 
state to the target state via TQD is tf = 72 /g when = 0.2 g, A = 2.3g, t 0 = 0.14 tj, 
t c = 0.19 tf and g = v, while the evolution time for the adiabatic passage is tf = 400 /g when 
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Q 0 — 0.2 g, t 0 = 0.14 tf, t c = 0.19 tf and g = v. So, the benefit of the TQD method is shown 
obviously that the speed via TQD method is faster than that via adiabatic passage. It is 
more worthy to note that the fidelity of the target state via TQD is almost equal to that 
via adiabatic passage. So our scheme has a huge advantage compared with the proposals 
via adiabatic passage. That means the present scheme via STAP method is not only fast 
but also robust. 

As we all know, it is necessary for a good scheme to tolerate the deviations of the ex¬ 
perimental parameters, because it is impossible to avoid the operational imperfection in 
experiment. Define that Sx = x' — x is the deviation of the ideal value x, x' is the actual 
value. In Fig. 10, we plot the fidelity of the target state | ^target) versus the deviations of the 
experimental parameters g, v, Do, and T (T — tf denotes the operation time). Numerical 
results demonstrate that our scheme is robust against the fluctuation of the experimental 
parameters. 


V. THE GENERATION OF THE A-ATOM GHZ STATE VIA 
TRANSITIONLESS QUANTUM DRIVING 

Next we briefly present the generalization of the scheme in Section IV to generate N- 
atom GHZ state by the same principle. We consider the physical configuration shown in 
Fig. 11, where N atoms cq, a 2 , ••• , cqv are trapped in N cavities Ci, C' 2 , ••• , Cn 
connected by N — 1 fibers / 1 , / 2 , • • • , /w-i, respectively. The level configurations of the 
atoms between two ends are the same as that of the atom a 2 in the three-atom case, and 
the level configurations of cq and a n are the same as those of cq and 03 in the three-atom 
case, respectively. The Hamiltonian of the present system can be written as in the rotation 
framework 


Htotal = H[ + H ', 

H'i = ^t'i\e) ai (g 0 \ + Q' N \ e )a N (9o\ + H.c., 

N -1 N 

H c ^ ^ \^)at {dl I T ^ ] 9j,r^j,r \ a, {dr 

i= 1 3 =2 
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N-l 

+ ['Ufc-i^_ 1 (ofc-i,/ + cik,i ) + Vkb\{ak, r + Ofc+i,r)] + H.c.. (28) 

k=2 

Let us consider the situation where N is an odd number, i.e., N — 21 + 1, (/ = 1, 2, 3, • • •). 
Suppose that the initial state of the atoms is \gogig r gig r • • • 5v) while all the cavities and 
fibers are vacuum, then the system can be expended in the following subspace 


101) = \gogwr ■ 

<7r) |0) a ZZ, |0 2 ) \&9l9r 

fl , r)|0)a/Z) 

|0' 3 ) = \gi9Wr ■ ■ 

•• 9r ) 11)cn |0l) = \9l9l9r •• 

flv)|l)/i> 

|0' 5 ) = \9l9l9r ■ ■ 

•• ^r)|10) C2 , |0' 6 ) = •• 

' 9r) 0) a ZZ, 

Wl) = 1 9Wrgr ■ 

<7r) 101) C2 , 10g) 1 9l9r9r 

fi'r) |l)/ 2 : 

|0' 9 ) = |0t0r0r • 

" 9r) 101)c 3 , 101 O ) = 1 9i9re ■ 

' ' <7r) |0) aZZ 


|0lv-i) = \9i9r ••• 0r0o)|O)aZI, (29) 


where |0) a ;; means that there is none photon in all boson modes, |nin 2 ) Si (s — C, f-i — 
1, 2, • • • , N) means that there are ri\ left-circularly photon and n 2 right-circularly photon 
in the corresponding cavity (7* or hber f % . 

Similar to the above procedure from Eq. (10) to Eq. (16), we get an effective Hamiltonian 
Heff(N) = KWMm + + H.C .), (30) 


where 


N 


N-l 


i« = iA- 2 > - E fw<». 


(31) 


i =1 2—1 

In addition, the eigenstates and eigenvalues of the Hamiltonian in Eq. (30) can be written 
as 


\Xo{t)) = 


( cos 6'{t) ^ 
0 

^-sin 9\t) ) 


, l.v.(O) = 


y/2 


( sin 6'(t) ^ 

±1 

l cos 0'(t) J 


(32) 
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with the corresponding eigenvalues x'o = 0 and X± — ±-V(fT, where tan O' = and 
hi' = + £1%. Substituting Eq. (32) in Eq. (19), we obtain 

H'(t)=id'\<i>' 1 ){<j>' iN -i\+H.c., (33) 


where O' = [fi' 1 (t)fl( v (t) — f2 / 1 (t)0( v (f)]/f2' 2 . 

Inspired by the above idea in section IV, we make the system into a non-resonant system 
to construct the CDD Hamiltonian in Eq. (33). Therefore, the Hamiltonian of the present 
system reads H' total = H[ + H' c + H' d , where H' d = YliLi A|e)(e|. Similar to the above 
procedure from Eq. (1221) to Eq. (1231) in Section IV, we obtain the final effective Hamiltonian 


H\ 


Q, 


■/2 


O' 


fe(N) 




Q r Q 


N 


3A 


3A 


3A 


(l < / ) l)(^ > 4Af-ll + l^4iV—11) (^11) ■ (34) 


For simplicity, we set 0 1 = fl N = O , the front two terms caused by Stark shift can be 
omitted and the Hamiltonian becomes 


H N = tt x (tM)Wur-i\+H.c., 


(35) 


where 0(,(t) = — To guarantee their consistency, we put a change that —> —i^N- 

Then the eigenstates of H e ff( n) become 


i m) 


( cos 8'(t) ^ 
0 

^ i sin 0'(t) y 


I x' ± (t)) 


1 

71 


sin0'(t) ^ 

±1 


V 


—i cos 0'{t ) y 


(36) 


and the CDD Hamiltonian H{t ) becomes 


H\t) = -efW 1 w AN _ 1 \-ef\4f AN _ 1 m 


(37) 


Compared Eq. (35) with Eq. (37), we can easily get the CDD Hamiltonian when the 
condition = —O' is satisfied. 


n[(t) = n' N (t) = n\t) = vTa^. 


(38) 


VI. EXPERIMENTAL FEASIBILITY AND CONCLUSIONS 

Now experimental feasibility needs to be discussed. The configuration of 8 ' Rb can be 
suitable for our proposals. Under current experimental condition a set of CQED parameters 
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g = 2n x 750 MHz, 7 = 2n x 2.62 MHz, and k c = 2n x 3.5 MHz are available with the 


wavelength in the region 630 — 850nm 


By using fiber-taper coupling to high-Q silica 


microspheres the efficiency of fiber-cavity coupling is higher than 99.9% 39]. The optical 
fiber decay at a 852nm wavelength is about 2.2dB/km 40], which means the fiber decay 
rate is about Kf — 1.52 x 10 5 i7A. With the above parameters, we obtain a relatively high 
fidelity about 97.15%. 

In conclusion, we have proposed an efficient scheme to fast deterministically generate 
TV-atom GHZ state in separate coupled cavities via transitionless quantum driving (TQD) 
only by one-step manipulation. We apply a promising method to construct STAP by joint 
utilization of the Zeno dynamics and the approach of TQD in the cavities QED system. 
The method features are that we do not need to control the time exactly and the evolution 
process is fast. Because the atoms are trapped in separate coupled cavity, the single qubit 
manipulation can be realized easily. When considering dissipation, we can see that the 
method is robust against the decoherences caused by the atomic spontaneous emission and 
fiber decay. The results show that the scheme has a high fidelity and may be possible to 
implement with the current experimental technology. So, the scheme is fast, robust and 
effective. We hope the scheme can be used to generate multi-atom GHZ state in the future. 
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FIG. 1: The structure of the experimental setup and atoms. Three identical atoms a\, 0,2 and 
are trapped in three separated cavities C\, C 2 and C 3 , which are linked by two fibers / 1 , f 2 - 



FIG. 2: The laser pulses and ^ 3/^0 versus t/tf. 
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FIG. 3: The structure of the experimental setup and atoms for the APF Hamiltonian. 



FIG. 4: The fidelity versus the laser pulses amplitude Hq and the operation time t/tf. 
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(a) (b) 


FIG. 5: (a) The population Ptarget of the target state | ^target) and the population Po of the 
initial state ^(O)) governed by the original Hamiltonian Hf 0 t via the adiabatic passage, (b) The 
population P\{t) of the states \4>\) and the population P\ \ (f) of the states governed by the 
original Hamiltonian Pitot via adiabatic passage. The parameters are collectively with the fixed 
values Hq = 0.2 g, g = v, to = 0.14 tf, t c = 0.19tf and tf = 400/#. 



FIG. 6: The fidelity F of the target state \ip(tf)) governed by Pl' tot versus the interaction time gtf 
and the detuning A /g. 
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FIG. 7: The population P\ (f) of the state \4>i) and the population P\ i (f) of the state |(/>n) governed 
by (a) the APF Hamiltonian H' tot with A = 2.3 g. (b) The original Hamiltonian Htot collectively 
with the fixed values Hq = 0.2 g, g = v, to = 0.14 tf, t c = 0.19 tf, and tf = 72/ g. 




(a) (b) 

FIG. 8 : The fidelity of the target state | i/j(tf)) governed by (a) the APF Hamiltonian H' tot with 
A = 2.3 g, tf = 72 jg and flo = 0.2 g. (b) the original Hamiltonian Htot with tf = 153 /g, and 
Oo = 0 . 5<7 collectively with the fixed values to = 0.14 tf, and t c = 0.19 tf versus the dimensionless 
parameters 7 / 5 , K c /g, and Kf/g, respectively. 
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FIG. 9: The fidelity of the target state | i/j(tf)} governed by (a) the APF Hamiltonian H' tot with 
A = 2.3 g, tf = 72/ g, and Ho = 0 fig. (b) the original Hamiltonian Htot with tf = 153 /g, and 
Ho = OAg 1 collectively with the fixed values to = 0.14tf, and t c = 0.19t/ versus the dimensionless 
parameters 7 /g and k c / g. 



FIG. 10: The fidelity of the target state | ^target) versus the deviations of (a) g and v, (b) T and 

Hq- 



tively trapped in IV-cavities which are linked by N — 1 fibers. 














